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1 Introduction 



Thermal Field Theories (TFT) represent a very wide subject and only a few aspects will be 
touched upon in these notes. The physical systems considered will be the hot QCD (quark- 
gluon) plasma or the hot QED (electron-photon) plasma in equilibrium. To probe the prop- 
erties of the plasma one often studies rare processes in the plasma which do not disturb the 
thermal equilibrium. As an illustration of such processes we shall consider the production of 
real or virtual (lepton pair) photons in a QCD plasma. The rate of production is proportionnal 
to e^, where e, the electric charge, is taken to be e <C g, with g the strong interaction coupling 
which is characteristic of forces maintaining the plasma in equilibrium. 

The following discussion will be based on perturbation theory and its improvements. The 
main problem one encounters in such an approach is that of divergences, both of the infra- 
red and the coUinear types, related to the masslessness of quarks and gluons. First the 
Feynman rules at finite temperature will be introduced and several versions of them will 
be presented depending on the formalism used: imaginary-time formalism (ITF), real-time 
formalism (RTF), retarded/advanced (R/A) formalism. The production of a virtual photon 
in a QCD plasma will be discussed and contrasted with the case at zero temperature. The 
improvement of perturbation theory when soft momenta are involved is presented: this is the 
hard thermal loop (HTL) resummation of Braaten-Pisarski and Frenkel and Taylor and some 
consequences are discussed. Further problems associated to mass singularities are discussed 
and the need to go beyond the HTL scheme is stressed. 



2 Feynman rules 

Thermal field theory is the application of the technics of usual quantum field theory which 
describes the interactions among a few fundamental "particles" , to the study of systems char- 
acterized by a large number of particles (statistical systems) at a given temperature (thermal 
equilibrium is assumed) [|I], 0, |^. The thermal expectation value of an operator A is defined 
by: 

<A>i3 = Z-^Tiie-f^^A) 

= Z^^^e-^-^™ < m|A|m > (1) 

m 

where the partition function is Z = Tr(e~^^), with f3 = 1/T the inverse of the temperature. 
If one remembers the time evolution for an operator in the Heisenberg representation 

e*^*'A(t) e-*^*' = A{t + t') (2) 

then one can write 

e^^A(t) e^^^ = A{t - t/3) (3) 

and interpret the inverse temperature as an imaginary time. In TFT it is then natural to intro- 
duce a complex time variable, the imaginary part of which being related to the temperature. 
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Consider now the thermal expectation value of a bilocal operator 
< A{t)B{t') >p = Z-^Ti (e-'^^A(t)5(t')) 

= <B{t' -zP)A{t)>^ (4) 

where the property of the cyclicity of the trace has been used. This equation summarizes the 
important Kubo-Martin-Schwinger (KMS) condition which expresses thermal equilibrium. 
Projecting on complete sets of states and introducing the evolution operators, the correlation 
function can be written as 

< A{t)B{t') >p = Z-^ J2 e'"'''" < m\A{t)B{t') \m > 

m 

= ^-1 ^ e-'^"(*-*')e'^'"(*-*'+'^) < m\A{0)\n >< n\B{0)\m > . (5) 

m,n 

For this expression to be defined, the exponentials should be well behaved when Em oo 
and this requires 

-P< Im(t - t') < 0. (6) 
We turn now to a scalar field theory where (f){x) is a real field. The propagator is 

< G,{x - x') >/3 = < 9,{t - t')(P{x)(P{x') + 9,{t' - t)^{x')4>{x) >p 

= <e,{t-t') G+{x-x') + e,{t' -t) G~{x-x')>p . (7) 

Remembering that the time variable can be complex, the time ordering operator T^, implied 
by the function 9c{t — t'), generalizes on an oriented contour C in the complex time plane, the 
time ordering operator defined on the real axis. Then, according to eq. (|^) the oriented contour 
has to be descending. The application of the KMS condition to the thermal propagator leads 
to the relations 

G-{t-t') = G+{t-if3-t') and G+{t - t') = G- {t - t' + i(3). (8) 
One obtains the propagator and other Green functions from the generating functional 

Z[j] = Tr (^e-^-^T, exp (^t ^ d'^x3{x)^{x)^^ (9) 

by the usual differentiation formula 

G.(x,......„) = -L^q£(£|L (10) 

Z[0\ tdj{xi) ■ ■■idj{Xn) 
Introducing the eigenvalue of the field operator at some initial time ti 

, ti) >= 0(f) ti) > (11) 
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one has to evaluate 

Z[j] = Jd(f)< 0,ti|e-^^T,exp (^t Jj^xj{x)(j){x)j > 

= J d(j)< (f),ti-ip\Tcexp(^i Jj^xj{x)^{x)^\(l),ti> . (12) 

Using well-known technics of splitting the time interval in small intervals one re-writes the 
generating functional as a functional integral 

Z[j] = N J [V{(j))]expi Jj'^x {C{x) + j{x)(t){x)) (13) 

where the time contour runs from tj to tj — i/3 and the fields are subject to the periodic 
boundary condition 

(j){ti,x) = (f){ti-i(3,x) (14) 

because of the same "initial" and "final" state as implied by the Trace operation. From now on, 
one constructs the perturbative expansion, as in the zero temperature case, by decomposing 
the lagrangian into its free part and its interacting part. 

We concentrate now on the construction of the propagator before specifying the contour 
which will lead to the different formalisms mentioned in the introduction. First, introduce the 
Fourier transform of the functions 

D^{K) = J d^xe'^'^'Cfix). (15) 

The KMS conditions easily give the relation 

D-{K) = e^'">D+{K) (16) 
and defining now the spectral function p{K) = D^{K) — D^{K) one derives 

D-{K) = p{K)n^{ko), withn^{ko) ^ 



e/3fco _ 1 

D^iK) = p{K){n,{ko) + l) (17) 



and in the space-time representation 



Gc{x) = J ^e-%(ir)(^,(t) + n,{ko)). (18) 

To obtain the free field propagator one expands the field into the creation and annihilation 
operators: 

0(x) = / ^ 2n9{ko)S{K' - m') [a(ir)e-^^^ + aHK)e'^^] (19) 

and one constructs the propagator from its definition eq. (|^). Making use of the usual 
commutation relations, e. g., 

[a{K), a\K')] = {271^ 2uJk 6{k - k') (20) 
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one finds the general expression 

^^^""-y^ = J ^se-^^^^-yh{k')6{K'-m'm{x'-y')+n,{ko)). (21) 

One immediately obtains the spectral function of the free propagator to be 

p{K) = 27Te{k'')6{K^ -m^). (22) 

To proceed further with actual calculations it is necessary to choose a contour and this will 
lead to the different formalisms with different expressions for the Green's functions. Of course, 
the predictions for rates or cross sections should be independent of the choice of the contour. 



2.1 The imaginary- time formalism 

This is based on the simplest contour, running vertically in the time interval [0, —if3] Jl], @, §]. 
It is convenient to introduce the "imaginary" time r, related to the "real" time hj t = —it. 
Using the KMS conditions one shows that the propagator is a periodic function of the energy of 
period 27r«//3. The propagator is expanded as a Fourier series with the Matsubara frequencies 
ujn = 27m/ P and one has 



or equivalently 



G{r,x) = -Y: /T^e— "V'=-A(^u;„,f) (23) 



A(2cj„,f) = dT j rfW""e-*'=^"G'(r,x). (24) 

^0 ^,0 







Using eq. (^) with 6c{x — y ) = 1 when < t < (3 one finds: 



[i) propagator: 



A^iun, k) = — — — 2 "with u;^ = P + m^; 



(ii) vertex: —ig as at T = 0; 

(iii) energy-momentum conservation: {2tiY5{K) — > — i/55„_o(2vr)^5(/c); 

(iv) loop integral: 

d^K i r d^k 

When working in the ITF an analytical continuation to real external energy variables must 
be performed 

iujn ^ k^ ±ie, e ^ 0+, (25) 
corresponding to "retarded" and "advanced" energies. 
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2.2 The real-time formalism 



The idea is to include the real axis in the oriented time contour as shown in Fig. 1 0, Q . The 
arbitrary parameter a is such that < a < /5. It is convenient to introduce the following 
propagators depending of the relative position of x and y (see Fig. 1 for the label of the various 
pieces of the contour): 



(i) 


if 




y on 


Ci 




G{x 


-y) 


= G'' 




-y) 


(ii) 


if 




y on 


C2 




G{x 


-y) 






-y) 


iii) 


if 
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on Ci 


, y on 


C2 


G{x 


-y) 


= G'^ 




-y) 


iv) 


if 


X 


on C2 


, y on 


Ci 


G{x 


-y) 


= G2i 




-y) 



G {x^ — y^ + i(T,x — y) 

G^i^x'^ — ia — ,y — x) 

The cases when one (or both) of the time coordinates is on a vertical part of the contour 
are not explicitely needed when working in momentum space although the vertical pieces of 
the contour play a crucial role in deriving a consistent formalism |^. From the general 
expressions above one obtains the Fourier transforms: 

D'\K) = - \ + 2vrn,(|A:°|) 6{K' - m') = (D'^k))* 

K'^ — — le 

D^\K) = 27r(^(-A;°) + n^(|A;°|) 5{K^ - m') e"^" 

D^\K) = 27r(0(A;°)+n^(|A;°|) (5(K2-m2) 6-'^'=°. (26) 

Similar formulae can be derived for a fermion and they involve the Fermi-Dirac fonction 
np{\k^\) = l/(e^'^°' + 1) rather than n^(|/c°|). A nice feature of the RTF is that it separates 
the temperature dependence contained in the statistical factors from the T = part. In the 
RTF there is effectively a "doubling" of fields since one may introduce 

(f)^{x) = (f)(x) when x is on the horizontal contour Gi 

4>^{x) = (f){x — ia) when x is on the horizontal contour C2. (27) 

This leads to two types of vertices: fields of type 0^ couple together with strength —ig or 
fields of type 0^ couple with strength ig, the change of sign in the latter case arising because 
the time variable is running from right to left. The propagator in the RTF is a 2 x 2 matrix 



^{k) = (D^^(k)) with elements defined in eqs. (pBf). When calculating Green's functions in 
this formalism it is necessary to sum over "type 1" and "type 2" internal vertices. This is 
crucial to cancel mathematically undefined terms, such as products of 5-functions with the 
same argument, which appear in intermediate stages of the calculation: to illustrate this point 
it is sufficient to consider self-energy insertions on a propagator and check that, only after 
summation over all types of internal vertices, terms of the form (5(7^^ _ 77^2) )"■ disappear. 

"Cutting rules" for calculating the imaginary part of Green's functions have been derived 
which are extremely useful to obtain rates or cross sections §]. The relation between 
RTF Green's functions and the analytic continuation of ITF ones is neither trivial nor simple 



[T0|, 0, |T2[ and has led to some confusion in the past [|1^, |T^. A nice way to relate the two 
in some cases, is described next. 
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2.3 The R/A formalism 



This formalism is a real-time formalism but it allows to construct directly retarded/advanced 
Green's functions, i.e. functions with external energies {ko + ie) or {ko — ie) which are related 
to the ITF Green's functions analytically continued as in eq. ( pSj) [ p!5| , |1^. The basic step is 
the following decomposition of the RTF propagator 



where the retarded and advanced propagators are given by 



V{K) (2^ 



.-0. (29) 

The matrices U and V depend on momentum K as well as the thermal factor and the pa- 
rameter a. They are associated to the RTF vertices and after summation on the two types of 
vertices one is led to the following momentum dependent vertices 

g,,AP,Q,R)=9nnniP,Q,R) = 0, 

9rra{P^ Q, P) = 9arr{P^ Q, P) = 9rar{Pi Q, P) = 9, 
g^,,{P,Q,R) = -g (1 + n,(gO) + n,(r°)), 

<7,,,(P,Q,i?) = -(7(l+n,(/)+n,(r°)), ••• (30) 

where the indices refer to the retarded/advanced prescriptions for the incoming external mo- 
menta P, Q and R. The first two equalities above express the KMS conditions in the R/A 
formalism. In loop calculations the sum over R and A internal indices is implied. Some useful 
discontinuity formulae to calculate cross sections have been derived in this formalism [15, 17]. 



3 Lepton pair production in a QCD plasma in pertur- 
bation theory 

It is instructive to consider first the production of a lepton pair in a QCD plasma: this 
illustrates some of the features and problems occuring in thermal calculations [l^, |20[| . 
When calculating this rate, at the first order in QCD, in hadronic collisions {i.e. at T = 0) it 
is well known that, at the partonic level, the infrared divergences (soft gluon) cancel according 
to the Lee-Kinoshita-Nauenberg theorem while the coUinear divergences survive due to the 
collinear emission of a gluon by the annihilating quarks (process qq — > gj*) or the collinear 
splitting of an initial gluon in the process gq — > q'j*. However, the partons are confined 
and these divergences are absorbed, via the factorization theorem, in a re-definition of the 
structure functions of the hadrons which become scale dependent: in other words confinement 
shields these singularities. Since, in a quark-gluon plasma, the partons are deconfined, it is 
legitimate to ask what happens to the collinear singularities. 
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For a pair of invariant mass squared Q^, produced at rest, the rate is ^1 

dN 



a 



Imn%(Q). 



(31) 



where Imll^ fiiQ) is the imaginary part of the trace of the retarded photon self energy. It has 
been calculated at the first order in QCD in the RTF and the result was found to be finite, of 



the form |22, 23 



a 



Imn%(g) = -Njl-2n 



IT 



1 H 

TT 



5 + F(«' 

4 ^T' 



(32) 



In the above equation we have considered only one species of massless quarks of charge e; 
is the fine structure constant of the strong interactions and and are the usual colour 
factors. Except for an overall normalisation factor, the thermal dependence is contained in 
the function F. The calculation shows that both infrared and collinear divergences cancel in 
the T = as well as in the thermal piece ||18|, |19], ^ . This occurs for two reasons: the phase 



space available at finite temperature is larger than at temperature (one has to integrate over 
initial particle momenta) and processes such as qqg — > 7* are included together with the usual 
annihilation and Compton diagrams of the Drell-Yan process (in this respect the factor 3/4 
in the correction term of eq. (^) is exactly the factor appearing in the T = annihilation 
of e'^e~ — >■ qqg); secondly there exist relations among the Bose and Fermi statistical factors 
which express thermal equilibrium. For example, in equilibrium the rate of transition for a 
gluon of energy Ei to be absorbed by a fermion of energy E2 is the same as the rate for a 
fermion of energy Ei + E2 to decay into a gluon of energy Ei and a fermion (detailed balance 
principle). This leads to relations of type 

n^{E,)n^{E2)=n^{Ei + E2){l + n^{E,)-n^{E2)) (33) 

and other similar relations. The cancellation of divergences has been explicitely shown to hold 
in scalar theories up to 3 loops p4[ but no general proof exists (see however |25|). Recently 
it has been shown that the cancellation of singularities can also occur out of equilibrium ||2r 
The thermal correction factor has been calculated and it behaves as 

T2 j2 



TT I 



rp2 rp2 



as 



as 



> 
II 



00. 



(34) 



This result can be understood as follows: the factor reflects the size of the thermal phase 
space while the is necessary for dimensional reasons. When the lepton pair invariant 

mass becomes of order gT the perturbation series appears not to be well behaved as the 
correction term becomes as large as the lowest order one and when — > a divergence 
appears again. This illustrates a general feature of thermal field theories where perturbation 
theory breaks down when "soft" scales, i.e. scales of order gT, are involved. This is put 
in a rigorous way in the "Hard Thermal Loop" (HTL) resummation scheme of Braaten and 
Pisarski ^ and Frenkel and Taylor ^ to which we turn next. In this approach, it 
is natural to distinguish between "hard" scales, of 0{T), the typical energies of quarks and 
gluons in the QCD plasma, from the soft scales, of 0{gT), with the assumption g <^1. 
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4 Hard Thermal Loop Resummation 



Consider, as an example, the retarded self-energy of a quark of momentum P in a QCD 
plasma. At the one-loop order, it can be expressed as 



(35) 

If the external momentum P is soft, with components of 0{gT), this expression reduces to 



hI - Mpo + lo)) e(/ + /°)5((P + Lf) A^{L) 



S«(P) -'-^Jd..in^ (.) + (-)) / ^ ^ (36) 

where we have introduced the light-like vector L = (1,/) and u = \lo\. The dimensional and 
the angular part of the loop integration factorize and it comes out p7| : 

^niP)'^^lf^J—. (37) 

1 J l-K PL + le 

The factor results from the uo integration and it is 

< = o.-^- (38) 

The behavior arises entirely from the region where u ^ T, i.e. when the loop momentum is 
"hard", the "soft" contribution of the loop being suppressed by factors of g. When inserting 
the self-energy correction on the fermion propagator the effective propagator is obtained: 

It is clear from eq. ( ^7|) that when P is soft the self-energy correction is also 0{gT) and the 
modification implied by the loop is of the same order as the bare propagator. On the contrary 
when P ~ T the self energy insertion induces a correction of 0{g). As a consequence, for 
a consistent perturbative calculation, one should use the effective propagator eq. ( |39D when 
the momentum is soft while the bare one is sufficient when the momentum is hard. The 
pole in the propagator eq. ( p^ ) leads to complicated dispersion relations with two branches 
(quasi-particle excitations) and the fermion acquires an effective mass rriq at rest given by 
eq. (|38D ; for hard momentum, one of the dispersion relations |]3T| , ^] leads to a particle-like 
excitation of mass 2mq whereas the other branch (called plasmino) decouples exponentially. 
An important consequence of thermal corrections is that the self energy acquires an imaginary 
part in the space-like region as can be seen by evaluating explicitely 

2 

ImSfi(P) = j dLUiPL) (40) 
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If < then one immediatly finds ImE ~ gT wliile if > it is vanishing, as the (5-function 
has no support (an exact calculation would give g^T), and the quasi-particle interpretation is 
justified. The origin of the imaginary part when P^ < is the Landau damping effect which 
accounts for the scattering of quarks and gluons in the thermal bath. Such a mechanism has 
no equivalent at zero temperature where the imaginary part exists only above the threshold 
P2 > 0. 

In the same way one finds that the gluon propagator has a hard thermal loop contribution. 
In fact, in the Landau gauge, the effective gluon propagator takes the form |3^ (for 
simplicity, we do not specify the R/A indices) 

*D^'^(L) = -P^%L)A'{L) - P^%L)A\l) (41) 

where 

A^'^(L) = — (42) 

The tensors P^*^ (L) are the transverse and longitudinal projectors whose explicit expressions 
are found in The vanishing of the denominators implies different dispersion relations for 
the transverse and longitudinal polarisation states. One has the following limits: 

n,(^,r= 0) = nju;,r= 0) = = g^T^[N + iVf/2]/9 (43) 

which means that the gluon acquires an effective mass in the plasma. On the other hand, in 
the static limit, 

n^(cu = 0,r^O) = (44) 
n^(cj = 0,r^O) = 3mJ (45) 

implying static screening of the electric field but not of the magnetic field. 

Other Green's functions may receive hard loop corrections. For example, for the quark- 
quark-photon vertex (which we will need later) the HTL contribution is simply written as 



_„V,(P,Q.«)^-4/|4^ (40) 



where Q is the photon momentum and L = (1, /) is as in eq. (PB]). When the external momenta 
P, Q and R are soft the integrand is ~ l/{gTY and the one loop vertex is of 0{e) as the bare 
QED vertex. On the contrary, if at least one of the external lines is hard the effective vertex 
is suppressed by at least one power of g. One is naturally led to introduce an effective vertex 

f ,(P, Q, R) = -ze(7A + V^a(P, Q, R)) (47) 

to be used, instead of the bare vertex, whenever all P, Q and R are soft. 

Likewise it can be shown that the three-gluon vertex, and more generally the ra-gluon 
vertex receive hard thermal loop corrections, as do Green's functions with 2-fermions and 
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n — 2 gauge bosons. Furthermore, in the HTL approximation the hard thermal loops are 
gauge invariant and they satisfy QED-like Ward identities For example, for the 3-gluon 
and the q — q — g vertex, respectively. 



Q^*tx{P,Q,R) = -i{*S-\P)+* S-\R)) (48) 



In conclusion, when carrying out a perturbative calculation, to obtain a consistent result bare 
soft lines and vertices should be replaced by their effective counterparts: resummed propaga- 



tors defined in eqs. (|39| , |4T|) and effective vertices {e.g. eq. p7|)). In Feynman diagrams effective 
Green's functions will be indicated by a • to distinguish them from their bare counterparts. 
Two-loop corrections are not needed in this scheme since internal lines in self-energy and, more 
generally n— point functions, being hard, corrections to these will be suppressed by factors 
of g. These rules can be deduced systematically from an effective gauge-invariant lagrangian 



3^, |3^ . The application of this scheme has lead to tremendous progress in the calculation of 



observables and to the solution of long standing puzzles. 



5 Applications of the HTL resummation 



We start by considering again the production of a soft virtual photon at rest [^. As discussed 
before one should calculate the imaginary part of the vacuum polarisation diagram which, at 
lowest order, describes the annihilation of a soft quark- ant iquark pair in a plasma. Since all 
momenta in the problem are soft one should use effective propagators and vertices and the 
graph to be considered is shown on Fig. 2. There exist other loop diagrams but, in principle, 
they do not contribute when summing over the photon polarisation indices. In terms of 
scattering amplitudes the graph of Fig. 2 has a very rich structure. Taking the imaginary 
part implies cutting through the effective vertices and propagators. Correspondingly, the 
rate of virtual photon production involves convolution of the diagrams in Fig. 3 |3^: if both 



internal fermion lines are time-like (P^, R^ > 0) the photon is produced in the annihilation 
or the decay of 2 quasi-particles (pole-pole) interaction (a); if one of the fermion line is 
space-like the imaginary part of the propagator corresponds to Landau damping and give 
rise to the scattering process (pole-cut) in (b); if both fermion lines are space-like then we 
obtain the double-scattering process (cut-cut) of (c). Cutting through effective vertices lead 
to similar diagrams (interference of (c) and (d)). The effect of Landau damping is of uttermost 
importance and increases the rate of production by several orders of magnitude compared to 
the estimates of eq. (^) |]36| . 

Similarly the rate of production of a hard photon has been calculated and found to be finite 
in the HTL approach while it would be infinite at lowest order of perturbation theory |38[] : 



indeed, in the annihilation and Compton processes of Fig. 4a, the fermion propagator may 
diverge: for a static exchange, for example, one has to evaluate an integral of the form / dp/p 
which diverges for soft 3-momentum transfer p. In the HTL approach, the effective fermion 
propagator eq. ( |57D should be used leading to the integral / dp 8p/ (8p^+mq(4+7r^)) regularized 
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by the fermion thermal mass. The rate is then 



go 



dW 
dq^ 



27r2 



(49) 



The calculation is much simpler than in the soft photon case as no effective vertices is needed 
since the vertices are connected to two hard lines: in fact it is sufficient to calculate the 
diagram of Fig. 4b where one of the internal lines at least is hard since the photon momentum 
is hard. Many more applications of the resummation scheme to the calculation of physical 
observables have been done . 

There has been in the past much controversy concerning the gluon damping rate {i.e. the 
imaginary part of the gluon self-energy): a negative rate would mean undamped oscillations 
and a perturbatively unstable quark-gluon plasma. Until the HTL scheme was developped 
a wide range of gauge dependent results could be found in the literature, based on the one- 
loop approximation, with rates of different signs and strength These results were in fact 



incomplete since they used only bare soft propagators and vertices. Defining the damping 
rate of a transverse gluon at rest as 7^,(0) = Imlly (mg, 0)/2mg, and performing a consistent 
calculation in the effective theory involves the evaluation of the graphs of Fig. 5 with effective 
propagators and vertices The result comes out gauge invariant and positive: 



7t(0) = 7.(0) = 6.635 



247r 



T 



(50) 



Among other successes one should mention the 0{g ) contribution to eq. (^3]) E3] 



6 Problems in the HTL resummation 

The use of the effective theory of Braaten and Pisarski makes it possible to improve the infrared 
behaviour of the theory in a gauge invariant manner. However there remains a number of 
difficulties related to the lack of static screening of transverse gluon modes (eq. (0)) and/or 
the appearence of collinear singularities when external particles are on-shell or massless. Two 
examples will respectively illustrate this problem: the damping rate of a fast moving particle 
and the production of soft real photons. 

Consider a fast moving fermion with energy pq ^ T. Its damping rate is given in the effec- 
tive approach by the imaginary part of the diagram of Fig. 6 with the dominant contribution 
arising from a soft transverse gluon and a hard internal fermion. The result is a convolution 
of the transverse gluon spectral function {i.e. the imaginary part of A (L) in eq. (|4^) when 

< 0) and one finds 

j{P)r^g\Tf Idl f ^p,{lo,l)r^g\T f (51) 

Jsoft J —I Lq Jsoft L 

This infrared divergence arises from unshielded static gluon exchanges (see eq. (PD). It 
should be noted that this divergence appears because one evaluates the diagram for an on- 
shell fermion, otherwise the off-shellness of the external fermion would shield the infrared 



11 



singularity. In QCD this IR divergence can be cured by the introduction of a magnetic mass 
(of non-perturbative origin) of 0{g^T) a solution not applicable to QED where gauge 



invariance requires the magnetic mass to vanish. Several proposals have been considered to 
solve the problem. An interesting recent solution is based on the observation that higher loop 
diagrams with static transverse gluons contribute to the same order as the graph of Fig. 6. 
Such corrections have been resummed in QED, by calculating the fermion propagator in 



a static background field and, defining the damping rate as the inverse of the decay time of 
the propagator in space-time coordinates, one finds 

^^g'T\n{^). (52) 

Going back to momentum space the retarded fermion propagator *S^{pQ,p) appears to be an 
entire function with singularity at Impo ^ —oo. 

The origin of the mass singularity problem can already be illustrated on a very simple 
example Coming back to eq. (pTf ) and constructing Tr(7°Ilij(P)) it comes out 

Re Tr(7°S^(P)) ~ --mj In (^^^] (53) 

which diverges logarithmically for the massless on-shell condition. Such a singular behaviour 
occurs also in the photon or gluon polarisation tensor and higher order diagrams contribute 
to same order as the lowest order one [^. Divergences like those in eq. ( [53| ) appear when 



calculating the production rate of soft real photons |T^- The diagram to be calculated is 
that of Fig. 2 and taking the imaginary part implies "cutting" trough an effective vertex (see 
an example in Fig. 7) i.e. taking the imaginary part of eq. ( ^61 ) and therefore to evaluate an 
integral of type: 

, ml fdL f 1 6(PL) ml r dL ^ ; 6(PL) , 

where the constraint P = R — Q has been used in the last equation. For a real photon, 
Q is light-like as is L and therefore the above integral exhibits a coUinear singularity when 
P and R are coUinear since PL and QL then vanish at the same point. Therefore, when 
applying strictly the HTL approach and keeping bare propagators and vertices when hard 
momenta are involved the production rate of soft real photon in a QCD plasma appears not 
to be defined. Even more interesting is the bremsstrahlung contribution to this rate [|^, ^ 
which is supposed to vanish in the HTL approximation as it is contained in the gluon tadpole 
diagram: this diagram summarizes in fact the graphs of Fig. 8, where the internal momenta 
P and R are now hard and L is soft, and involves the evaluation of: 



Imn%(Q) ~ eV go / d^Pd'L (^^^ n'^{po)nMp,Jlo,l) ^^^^Ip^'iy^^ ^^^^ 

Using the 5-function contraints the term l/i?^(P + L)^ naively appears of 0{1/ g'^T^) for P, R 
hard and L soft. However the denominator R^{P + LY is responsible for coUinear divergences 
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which drastically modify this naive estimate. Using the constraints one easily rewrites 



R2(p + Ly ~ P-QP-Q + Q-L~ [p-Q^ P-Q + Q-l)"^ Q-LP-Q' ^^^^ 

The first equality shows the presence of two very close collinear singularities (when P ■ Q=0) 
since the two poles differ only by the soft Q ■ L term. The last equality holds true to leading 
order only after the integration over the whole phase space is performed. Introducing the 
angular variable u = 1 — cos 6 between the light-like momenta P and Q the above expression 
becomes, near u = 0, 

' ' (57) 



R'^{P + LY qL^ pqu 

This form shows the presence of a logarithmic collinear divergence and the order of the residue 
at the pole in u is instead of the naively expected The near overlap of the 

collinear singularites causes the enhancement of the predicted rate by a factor l/g^. The 
cross section is regularised by keeping the soft kinematic terms in the hard propagator of the 



diagrams as well as the hard fermion thermal mass |^ which enters the calculation at the same 
order. This discussion can be generalised to the case of the production of quasi-real photons 
at soft momenta and the relevant parameter is found to be Q'^/q^ [|^, 0: when it is small 
(< g^) the rate of production is enhanced by a factor l/g"^ compared to the expected order 
of magnitude in the HTL resummation program while when Q'^/q^ ~ 1 the bremsstrahlung 



contribution is of the same order as that estimated in 15^ 



The implications of light-cone singularities for the HTL resummation program present an 
interesting challenge and they are still to be investigated |ET . 



7 Acknowledgements 

I thank F. Gelis for many discussions and for much help in the preparation of these notes. I 
also thank the organisers of WHEPP 4 for a very enjoyable meeting. 



References 

[1] A. Fetter, J. Walecka, Quantum Theory of Many Particle Systems, McGraw-Hill, New- 
York, 1971. 

[2] J.l. Kapusta, Finite- Temperature Field Theory, Cambridge University Press, 1989. 
[3] R. Mills, Propagators for Many-Particle Systems, Gordon and Breach Science Publishers, 
New- York, 1969. 

[4] N.P. Landsman, Ch.G. van Weert, Phys. Rep. 145, 141 (1987). 
[5] T.S. Evans, Phys. Rev. D 47, 4196 (1993). 
[6] F. Gelis, Z. Phys. C 70, 321 (1996). 

[7] R.L. Kobes, G.W. Semenov, Nucl. Phys. B 260, 714 (1985). 
[8] R.L. Kobes, G.W. Semenov, Nucl. Phys. B 272, 329 (1986). 



13 



[9] R.L. Kobes, Phys. Rev. D 42, 562 (1990). 
[10] R.L. Kobes, Phys. Rev. Lett. 67, 1384 (1991). 
[11] T.S. Evans, Nucl. Phys. B 374, 340 (1991). 
[12] F. Guerin, Nucl. Phys. B 432, 281 (1994). 
[13] R. Baier, B. Pire, D. Schiff, Z. Phys. C 51, 581 (1991). 
[14] H. Nakkagawa, A. Niegawa, H. Yokota, Phys. Lett. B 244, 63 (1990). 
[15] P. Aurenche, T. Becherrawy, Nud. Phys. B 379, 259 (1992). 
[16] M.A. van Eijck, Ch.G. van Weert, Phys. Lett. B 278, 305 (1992). 
[17] P. Aurenche, T. Becherrawy, E. Petitgirard, ENSLAPP-A-452/93, Fep-ph/9403320 
[18] J. Cleymans, L Dadic, Z. Phys. C 45, 57 (1989). 
[19] R. Baier, B. Pire, D. Schiff, Phys. Rev. D 38, 2814 (1988). 
[20] T. Altherr, P. Aurenche, T. Becherrawy, Nucl. Phys. B 315, 436 (1989). 
[21] C. Gale, J.I. Kapusta, Nucl. Phys. B 357, 65 (1991). 
[22] T. Altherr, P. Aurenche, Z. Phys. C 45, 99 (1989). 
[23] Y. Gabellini, T. Grandou, D. Poizat, Ann. of Phys. 202, 436 (1990). 
[24] T. Grandou, M. Le Bellac, D. Poizat, Nucl. Phys. B 358, 408 (1991). 
[25] A. Niegawa, Phys. Rev. Lett. 71, 3055 (1993). 
[26] M. Le Bellac, H. Mabilat, Universite de Nice preprint INLN-96/17. 
[27] E. Braaten, R.D. Pisarski, Nucl. Phys. B 337, 569 ( 1990). 
[28] E. Braaten, R.D. Pisarski, Nucl. Phys. B 339, 310 (1990). 
[29] J. Frenkel, J.C. Taylor, Nucl. Phys. B 334, 199 (1990). 
[30] J. Frenkel, J.C. Taylor, Nucl. Phys. B 374, 156 (1992). 
[31] V.V. Klimov, Sov. Phys. JETP 55, 199 (1982). 
[32] H.A. Weldon, Phys. Rev. D 26, 2789 (1982). 
[33] H.A. Weldon, Phys. Rev. D 26, 1394 (1982). 
[34] J.C. Taylor, S.M.H. Wong, Nuc. Phys. B346, 115 (1990). 
[35] E. Braaten, R.D. Pisarski, Phys. Rev. D 45, R1827 (1992). 
[36] E. Braaten, R.D. Pisarski, T.C. Yuan, Phys. Rev. Lett. 64, 2242 (1990). 
[37] J.I. Kapusta, P. Lichard, D. Seibert, Phys. Rev. D 44, 2774 (1991). 
[38] R. Baier, H. Nakkagawa, A. Niegawa, K. Redlich, Z. Phys. C 53, 433 (1992). 
[39] M. Le Bellac, Thermal Field Theory, Cambridge University Press, 1996. 
[40] R.D. Pisarski, Nucl. Phys. A 525, 175 (1991). 
[41] E. Braaten, R.D. Pisarski, Phys. Rev. D 42, 2156 (1990). 
[42] H. Schultz, Nuc. Phys. B. 413, 253 (1994). 
[43] A. Rebhan, Phys. Rev. D 46, 482 (1992). 
[44] R.D. Pisarski, Phys. Rev. D 47, 5589 (1993). 
[45] J. P. Blaizot, E. lancu, Phys. Rev. Lett. 76, 3080 (1996). 
[46] F. Flechsig, A.K. Rebhan, Nucl. Phys. B 464, 279 (1996). 
[47] R. Baier, S. Peigne, D. Schiff, Z. Phys. C 62, 337 (1994). 
[48] P. Aurenche, F. Gehs, R. Kobes, E. Petitgirard , |hep-ph/9606287 . 



[49] P. Aurenche, F. Gehs, R. Kobes, E. Petitgirard , |hep-ph/9609256 



[50] F. Gehs, preprint ENSLAPP-A-613-96, |IIep-ph/ 9608425 . 



[51] A. Rebhan, preprint TUW 96-15, |hep-ph/9608"4T0 



14 



Figure Captions 

Fig. 1 The real-time contour. 

Fig. 2 Diagram in the resummed theory contributing to the production of soft photon in a 
quark- gluon plasma. All momenta are soft. 

Fig. 3 The various processes involved in the evaluation of the diagram of the previous figure. 

Fig. 4 (a) Processes contributing to hard photon production in a quark gluon plasma; 
(b) The corresponding diagram in the Braaten-Pisarski effective theory. 

Fig. 5 Diagrams contributing to the gluon damping rate. 

Fig. 6 The dominant diagram contributing to the damping rate of a fast fermion; the gluon in 
the loop is transverse. 

Fig. 7 One of the diagram leading to a divergent contribution in the production of a soft real 
photon. 

Fig. 8 The dominant bremsstrahlung diagrams for the production of a soft real photon. 
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